Hydrodynamics with a chiral hadronic equation of state including quark degrees of 

freedom 
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We investigate the influence of a deconflnement phase transition on the dynamics of hot and dense 
nuclear matter. To this aim a hybrid model with an intermediate hydrodynamic stage for the hot 
and dense phase of the system is employed for collisions of Pb+Pb/Au+Au at beam energies of 
^^lab = 2 — IQQA GeV, while initial and flnal interactions are performed by a microscopic transport 
approach (UrQMD). In the hydrodynamic stage an equation of state that incorporates a critical 
end point (CEP) in line with lattice data is used. It follows from coupling the Polyakov loop (as an 
order parameter for deconflnement) to a chiral hadronic SU(3)f model. In this conflguration the EoS 
describes chiral restoration as well as the deconflnement phase transition. We compare the results 
from this new equation of state to results obtained, by applying a hadron resonance gas equation of 
state, focusing on bulk observables deemed to be sensitive to the phase transition to a Quark-Gluon 
Plasma. 
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Heavy ion collisions at intermediate incident beam 
energies (i?iab = 5 — 200^ GeV) offer the unique 
opportunity of being able to scan a wide range of 
temperatures T and baryo-chemical potentials /is in the 
phase diagram of strongly interacting matter 1] (for re- 
cent lattice QCD results see 0-11], for phenomenological 



studies see |5l-lll|). In this energy region one hopes to 



find experimental evidence for a deconfincmcnt phase 
transition from hadronic matter to the Quark-Gluon 
Plasma (QGP) phase (where quarks are deconfined). 
Especially the so-called critical end point (CEP), a 
point in the phase diagram that terminates the first 
order phase transition-line (which is expected for high 
chemical potentials), is of great interest. 



Key bulk observables like the directed flow vi 1^ 24 1. 
but also particle multiplicities, ratios, and their fluctua- 
tions, have been predicted and sometimes already shown 
to be sensitive to the active degrees of freedom in the 
early stage of the reaction. Indeed, the energy depen- 
dences of various observables show anomalies at low SPS 
energies which might be related to the onset of decon- 
flnement and chiral symmetry restoration . 

Early-on, fluid dynamics has been proposed as an el- 
egant way to include the EoS of strongly interacting 
matter in the description of heavy ion collisions (25l - 
27[. Especially since experiments at the RHIC facility 
have claimed to have found a (s)QGP that behaves like 
a nearly ideal fluid, the idea of modeling heavy ion col- 
lisions with (ideal) fluid dynamics has been revived 28- 

In order to study the fluid dynamical evolution of a 
heavy ion collision, the different boundary conditions (i.e. 
the initial space-time distributions of the corresponding 



energy and baryon density, as well as the freeze out pre- 
scription) have to be determined. Since experimental 
data provide mainly information from the final state of 
the reaction, integrated over the time evolution of the 
system, the initial state for hydrodynamical simulations 
is usually inferred from model assumptions or by an 'ed- 
ucated guess' in comparison to data. The connection be- 
tween (observed) final state and the inferred initial condi- 
tions is further blurred by the unknown equation of state, 
potential viscosity effects, and problems in the freeze-out 
treatment. Another issue concerns the assumption of 
thermal equilibrium, which is probably not fulfilled for 
the early stages of heavy ion collisions at intermediate 
energies. 

There have been attempts to solve these problems 
by describing such collisions with viscous or multi-fluid- 
hydrodynamic models 3ll-l37j|. but the practical applica- 
tion of these models is difficult. 

Transport theory offers another, different, approach 
aiming at the consistent description of heavy-ion reac- 
tions, from the initial state to the flnal decoupling of 
the system. This microscopic description has been ap- 
plied quite successfully to the partonic as well as to the 
hadronic stage of heavy ion collisions 38-4^. However, 
to explain hadronization and the phase transition be- 
tween the hadronic and the partonic phase on a micro- 
scopic level is one of the main issues to be resolved. It is 
therefore difficult to find an appropriate prescription of 
the phase transition in such a microscopic approach. 

To obtain a more comprehensive picture of the 
whole dynamics of heavy ion reactions various so called 
micro-|-macro hybrid approaches have been developed 
during the last years |41|. The NEXSpheRIO approach 
uses initial conditions that are calculated in a non 
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equilibrium model (NEXUS) followed by an ideal 
hydrodynamic evolution [42l444| . For the freeze-out a 
continuous emission scenario or a standard Cooper-Frye 
calculation is employed. Oth er g roups, e.g. Hirano et 



al 
or 



45l-l47t , Bass/Nonaka |48l - l50l |. use smooth Glauber 
CGC initial conditions followed by a full three- 
dimensional hydrodynamic evolution and calculate the 
freeze-out with a subsequent hadronic cascade. The 
separation of chemical and kinetic freeze-out and final 
state interactions like resonance decays and rescatterings 
are taken into account. 

In this paper, we apply a transport calculation with an 
embedded three-dimensional ideal relativistic one-fluid 
calculation for the hot and dense stage of the reaction, 
thus reducing the parameters for both the initial con- 
ditions and the freeze-out prescription. This will allow 
us to compare calculations with different EoS within the 
same framework and to extract the effect of changes in 
the EoS - e.g. a phase transition from hadronic mat- 
ter to the QGP - on observables. In [sH this model has 
been used, applying a hadron resonance gas EoS, to pro- 
vide a baseline calculation, disentangling the effects of 
the different assumptions for the underlying dynamics in 
a transport vs. hydrodynamic calculation. We will ex- 
tend this purely hadronic calculation by introducing an 
EoS that includes a deconfinement phase transition. The 
calculations will be performed in the broad energy range 
from E\ah — 2 — IGOA GeV where experimental data from 
BNL-AGS and CERN-SPS exists and which will be ex- 
plored in more detailed energy scans by the FAIR project 
near GSI and the RHIC low energy program. 



THE HYBRID MODEL 

The Ultra-relativistic Quantum Molecular Dynamics 
Model [52, l53| (in its cascade mode) is used to calculate 
the initial state of a heavy ion collision for the hydro- 
dynamical evolution [13]. This is done to account for 
the non-equilibrium dynamics in the very early stage of 
the collision. In this configuration the effect of event- 
by-event fluctuations of the initial state is naturally in- 
cluded. The coupling between the UrQMD initial state 
and the hydrodynamical evolution happens at a time 
^start when the two Lorcntz-contractcd nuclei have passed 
through each other: 
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(1) 



where R is the radius of the lead nucleus and 7c. m. 
the Lorenz gamma factor of the two colliding nuclei in 
their center of mass frame. At this start time all ini- 
tial collisions have proceeded, i.e. also the initial baryon 
currents have decoupled from each other, and it is the 



earliest time at which local thermodynamical equilibrium 
may be achieved. To map all 'point-like' particles from 
UrQMD onto the spatial grid of the hydrodynamic model 
each hadron is represented by a Gaussian of finite width. 
This procedure is necessary, since the cell length of 0.2 
fm is much smaller than the actual size of a hadron. The 
Gaussian width is chosen to be a = 1 fm. This width 
reflects the typical size of hadrons and avoids numer- 
ical instabilities (numerical entropy production) in the 
initial phase of the hydrodynamical evolution. This in- 
stantaneous thermalization at ^start goes along with an 
increase in entropy, as entropy is maximized in the equi- 
librium state. However, note that is has been checked 
that the results only weakly depend on the choice of the 
time tstart [EH . Especially final particle multiplicities and 
their mean transverse mass only change by about 10%, 
if the time tstart is doubled. 

For calculations at finite impact parameter, the spec- 
tators - particles that have not interacted until tstart - 
are propagated separately from the hydrodynamic evolu- 
tion. They are treated as free streaming particles until 
the end of the hydrodynamic phase has been reached. 

The full (3-1-1) dimensional ideal hydrodynamic evolu- 



tion is performed using the SHASTA algorithm [55|, [5 



The partial differential equations are solved on a three- 
dimensional spatial Eulerian grid with fixed position and 
size in the computational frame. The size of the grid is 
200 cells in each direction, while the cell size has been 
chosen to be = 0.2 fm which leads to time steps of 
dt = 0.08 fm in order to avoid non-causal effects in the 
propagation (Courant criterion). 

To transfer all particles back into the UrQMD model, 
an approximate iso-eigentime transition is chosen (see 
57| for details). Here, we 'freeze out' individual trans- 
verse slices, of thickness Az = 0.2fm, at a constant 
time-like transition hypersurface. This time for each 
slice is given, whenever the energy density e, in every cell 
of this slice, has dropped below four times the ground 
state energy density (i.e. ^ 580MeV/fm^). This assures 
that all cells have passed through the mixed phase of the 
equation of state and the effective degrees of freedom at 
the transition are hadronic (see Fig. [T] for a depiction 
of this line of constant energy density in the T-fig phase 
diagram) . 

By applying a gradual transition one obtains an almost 
rapidity independent switching temperature. The 
hydrodynamic fields, in a given slice, are transformed 
to particle degrees of freedom via the Cooper-Frye 
equation on an isochronous time-like hypersurface in 
the computational frame (the hypersurface normal is 
da^ ^ (d^x, 0,0,0)). 

As different longitudinal slices have different freeze 
out times, dcr^ should of course also have a space like 
component. Such a parametrization of the hypersurface 
is not easily dealt with numerically (especially since our 
system has locally fluctuating densities and therefore 
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FIG. 1: (color online) Isentropic expansion paths (red lines) 
in the T — ^q plane for very central Pb+Pb/Au+Au reactions. 
Isentropic expansion from the overlap model initial conditions 
are shown as full line in blue. Beam energies are from left to 
right: S^b = 40, 30, 20, 10, 5A GeV. The line of the first order 
phase transition is indicated in black together with the critical 
endpoint of the model. Also shown is the line of constant 
energy density e = 4eo (gray dashed) . 



an inhomogeneous hypersurface). We therefore try to 
justify our approach by comparing freeze out results 
from a simple, analytically solvable, one-dimensional 
Bjorken scenario by using both, our approximate and 
the correct parametrization of da^. In a realistic set 
up, where the longitudinal expansion of the system is 
about 10 fm, we obtain a total difference in particle 
production of about 10%. The error in particle pro- 
duction per rapidity interval grows for larger rapidities. 
In consequence, we expect the calculated rapidity 
distributions, of particles produced at the highest SPS 
energies (-Eiab — 160A GeV), to show the largest effect 
of our choice of dcr^j. More precisely, we expect the 
present rapidity distributions at the highest energies to 
be lower at mid-rapidity and broader at high rapidities, 
as compared to results with the correct dcr^. As a 
remark, it is possible to numerically extract the correct 
parametrization of the full hypersurface using digital 
image processing techniques [58|, and then compare 
these results with data acquired with our simplified 
hypersurface. There is work in progress on this task, 
which will be subject of future publications. 



As has been pointed out in Si], the present transition 
procedure conserves the baryon number, the electric 
charge and the total net strangeness on an event-by- 
event basis, but the total energy only on average and 
only if the freeze out is treated properly. In the present 
calculations, the total energy of the system varies by 
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FIG. 2; Binding energy per nucleon as a function of the net 
baryon number density at T = 0. The minimum at nuclear 
ground state density (po) corresponds to a binding energy of 
-16 MeV. 



1 — 2% per event. To overcome this problem it is possible 
to rescale the momenta of all produced particles in 
every single event to enforce exact energy conservation. 
We have employed both methods and checked that the 
results, for particle rapidity and momentum spectra, do 
not depend on the method used but are identical, within 
statistical errors, when averaged over 500 events. 

After the particles are created according to our pre- 
scription, they proceed in their evolution in the hadronic 
cascade (UrQMD) where rescatterings and final decays 
are calculated until all interactions cease and the system 
decouples. 

A more detailed description of the hybrid model in- 
cluding parameter tests and results can be found in fsij . 



THE EQUATION OF STATE 

After the introduction of the dynamical framework, 
we now turn to a detailed discussion of the novel EoS. 
Here, we employ a single model to obtain the EoS of the 
hadronic and the quark phase. 

The hadronic part of the model is a flavor-SU(3) model, 
which is an extension of a non-linear representation of a 
sigma-omega model including the lowest-lying multiplets 
of baryons and mesons (for the derivation and a detailed 
discussion of the hadronic part of the model Lagrangian 



see [59|-|6l[). In spirit similar to the PNJL model [62 1 
it includes the Polyakov loop $ as an effective field and 
it adds quark degrees of freedom. The temporal back- 
ground field $ is defined as $ = ^Tr[exp(i / drAi)], 
where A4 = lA^ is the temporal component of the SU(3) 
gauge field. 
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The Lagrangian density of the model in mean field ap- 
proximations reads: 



L — Lf^in + L. 



int I" -^mesonj 



(2) 



where besides the kinetic energy term for hadrons and 
quarks, the terms 



Lint = i-'ihoigiij^ + gi,i,(t>) + (3) 



'meson 2 ^ 

,4 I 0"^ I Q..,2a2 



-k2 (4 + C") - ha^C 
-ki In ^ + mlf^a 
+ (V2m2./, - ^^mlu) C , 



(4) 



represent the interactions between baryons (and 
quarks) and vector and scalar mesons the self interactions 
of scalar and vector mesons and an explicit chiral symme- 
try breaking term. The index i denotes the baryon octet 
and the three light quarks. Here, the mesonic conden- 
sates (determined in mean-field approximation) included 
are the vector-isoscalars uj and and the scalar-isoscalars 
a and C (strange quark- antiquark state). At this point 
we neglect the p-meson contributions as we only discuss 
isospin-symmetric matter. 

The effective masses of the baryons and quarks are 
generated by the scalar mesons except for a small explicit 
mass term and the term containing the Polyakov field $ 



(5) 



(6) 



With the increase of temperature/density, the scalar 
fields decrease in value, causing the effective masses of the 
particles to decrease towards chiral symmetry restora- 
tion. The Polyakov loop effectively suppresses baryons at 
high temperatures/densities and quarks at low temper- 
atures/densities due to their corresponding mass shifts 
shown above. 

Due to meson vector interactions, the baryons obtain an 
effective chemical potential: 



(7) 



All thermodynamic quantities are derived from the 
grand canonical potential 51 by assuming an equilibrated 
state, where —p = -p^ (p is the pressure) is in the absolute 
minimum of J7 with respect to all fields. For the region 



of phase coexistence (at the 1. order phase boundary) 
we assume a phase mixture. This of course implies, that 
the expanding system is in chemical equilibrium during 
the whole evolution. The grand canonical potential of 
the model has the form: 
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(8) 



Here Qth includes the heat bath of hadronic and quark 
quasiparticles (as function of T, ii* and to*) within the 
grand canonical potential of the system. The Polyakov- 
loop potential U will be discussed in the following. 
The potential U for the Polyakov loop reads: 



U = (aoT* ^ 
+a3ToMn(l 



8$3 - 3$^). 



(9) 



ft is based on [65|, |66| and fitted to the pressure and 
Polyakov loop values as computed in lattice-QCD calcu- 
lations at zero chemical potential as discussed in detail in 
67|. Additional terms, depending on the chemical poten- 
tial, are fixed in order to reproduce the phase diagram at 
high densities. This includes a first order phase transition 
line in yUg and T that ends in a critical point, of second 
order, at the values obtained by lattice calculations 



The coupling constants for the baryons (already shown in 
[6l[) are chosen to reproduce the vacuum masses of the 
baryons and mesons, nuclear saturation properties as well 
as the hyperon potentials. The vacuum expectation val- 
ues of the scalar mesons are constrained by reproducing 
the pion and kaon decay constants. The coupling con- 
stants for the quarks {gqa = —3.0, gs( — —3.0, Tq = 200 
MeV, ao = 1.85, ai = 1.44x10-3, 02 = 0.08, 03 = 0.40, 
gN<i> = 1500.00 MeV, gq^ = 500 MeV) are chosen to 
reproduce lattice data (for To = 270 and pure gauge a 
first order phase transition at = and T = 270 MeV 
is reproduced) and known information about the phase 
diagram. 

As can be seen in Fig. [T] the transition from hadronic 
to quark matter obtained is a crossover for small chemical 
potentials. At vanishing chemical potential the transition 
temperature is 171 MeV, determined as the peak of the 
change of the scalar field and the Polyakov loop. Beyond 
the critical end-point (at fi^^B = 354 MeV, 3^= 167 MeV 
for symmetric matter in accordance with [68|) a first or- 
der transition line begins. As can be seen in Fig. [2l the 
model reproduces nuclear matter saturation at realistic 
values for the saturation density, nuclear binding energy, 
as well as compressibility and asymmetry energy. In ad- 
dition, realistic results at low densities for the nuclear 
matter liquid-gas phase transition are obtained and the 
model has been successfully applied in order to model the 
properties of compact stars |63j . It is crucial for the nu- 
merical studies of heavy-ion simulations to have an equa- 
tion of state at hand that shows reasonable behavior over 
a large range of densities and temperatures, which is the 
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FIG. 3: (color online) Isentropic expansion paths (black lines) 
and contours of the speed of sound in the e-n phase diagram. 
The beam energies, associated with the lines, are the same as 
in the previous figure. The gray region relates to unphysical 
combinations of e and n (T < 0). The critical endpoint of the 
model is displayed as the black dot. 



FIG. 4: (color online) Excitation function of the averaged 
speed of sound in most central A+A collisions. The decon- 
finement EoS (red soUd line) is compared to the hadron res- 
onance gas (blue dashed line). Indicated is also the speed 
of sound extracted from 3 particle correlation studies at the 
STAR experiment [H]. 



case here. In the following we will refer to this EoS as 
the deconfinenient EoS (DE). 

We will compare results obtained with the DE to cal- 
culations using an EoS consisting of a hadron resonance 
gas including all reliably known resonances with masses 
up to 2 GeV (referred to as hadron gas HG). The HG 
is a very important ingredient of the model, because the 
active degrees of freedom on both sides of the transition 
hypersurface have to be equivalent to ensure the conser- 
vation of important quantities (e.g. entropy). For the HG 
this is the case, as it has the same degrees of freedom as 
the UrQMD model. In the DE the hadrons acquire effec- 
tive masses due to interaction and therefore this equiva- 
lence condition is only approximately fulfilled. To solve 
this problem we change the active equation of state after 
the last step of the hydrodynamical evolution (from the 
deconfinement EoS to the HG), thus obtaining the cor- 
rect temperatures and chemical potentials for the parti- 
cle distributions. If we use the DE EoS for the freeze-out 
prescription we find that the total energy conservation is 
violated systematically by about 3%. Therefore results 
on particle spectra also change slightly. As this change 
is only on the level of a few percent, we prefer to have 
the correct degrees of freedom during the transition from 
hydrodynamics to transport. 

Fig. [T] depicts the phase structure of the obtained 
EoS. Included are lines of constant entropy per baryon 
as they are expected for beam energies of E^iab — 
40, 30, 20, 10, 5A GeV (from left to right). The values for 
^ were calculated using a simple overlap model where: 
rib = 27c.m.'T-o, e = Vsjc.m.no, and Uq is the nuclear 



ground state baryon number density. It is obvious that 
the deconfined phase is already reached at energies above 
E\ah = 10v4 GeV. As has been shown in [5J], an incident 
beam energy of E\a,h — 40A GeV may not be sufficient 
to reach the critical endpoint suggested in [g^. If in- 
stead the CEP is situated at higher /ib > 450 MeV (as 
suggested by recent lattice studies [87|) it could again be 
accessible in the energy range of the FAIR project. 



THE SPEED OF SOUND 

An important property of a hot and dense nuclear 
medium is the speed of sound (c^): 



dp 



S/A 



dp 
de 



dp 
p dn 



(10) 



It is not only closely related to expansion dynamics but 
also controls the way perturbations (sound waves) travel 
through the fireball. Fig. |3] shows Cg as a contour plot in 
the e — n phase diagram, where e is the energy density 
and n the net baryon number density, both given in units 
of the ground state values (eo — 146 MeV/fm'^, no — 
0.159 fm-3). 

One can clearly see the reduction of the speed of 
sound (softening) of the EoS in the mixed phase. Also 
depicted is the location of the CEP and isentropic paths 
for energies as in Fig. [T] Please note that the numerical 
accuracy for Fig. 2 is very limited when going to very 
low temperatures. This is essentially because the energy 
density is proportional to T**. The T = line can only be 
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FIG. 5: (color online) Total Att multiplicities of pions (lower 
panel), H's (middle panel) and kaons (upper panel). Data 
[69h79i] are indicated by squares. 



thought of as a guide for the eye (the error in T is about 
30 MeV for the lowest value of e/eo). In fact we have 
checked that the binding energy and compressibility is 
indeed well reproduced at nuclear saturation density 
67| . In order to show that the system reproduces the 
correct nuclear groundstate, Fig. [2] depicts the binding 
energy per nucleon at T = 0, beeing of the order of a 
few MeV at p — pQ. 



To quantify this softening, we can calculate the aver- 
age speed of sound during the hydrodynamic evolution. 
We define the average {cs{t)) at a given time t as the 
average speed of sound over all fiuid cells weighted with 
the energy density of that cell, where the Cs of every cell 
can be deduced from the EoS as a function of energy 
and baryon number density. The speed of sound {cs(t)) 
is then averaged over the whole time evolution, where 
every time step has the same statistical weight. 

Fig. m shows the excitation function of the averaged 
speed of sound for both equations of state considered. 
The hadron gas (blue dashed line) yields higher values of 
(cs) than the deconfinement EoS (red solid line). This 
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FIG. 6: (color online) Rapidity distributions for different par- 
ticle species at three different beam energies. Hybrid model 
calculations, DE (solid line) and HG (dotted lines) results are 
compared to data (symbols) [t^. IssI]. 



is expected, since the phase transition leads to a soften- 
ing of the EoS. Still, in both cases the averaged speed of 
sound is well above 0.3 c. It was proposed, that a coni- 
cal Mach wave created by in medium jets traversing the 
hot and dense system of a relativistic nuclear collision, 
could provide the means to experimentally measure the 
speed of sound in the fireball. Indeed, experiments at the 
RHIC claim to have observed conical emission in heavy- 
ion collisions. Applying a 3-particle correlation method. 



the Mach angle 9m was extracted from data [GJ]. In a 
simple Mach cone picture this angle can easily be related 
to the speed of sound: 



cos (6* A/) = Cs/Vp, 



(11) 



where Vp is the velocity of the projectile creating the wave 
(usually Vp is considered to be close to the speed of light). 
This simple approximation leads to an estimate for the 
speed of sound of (c^) ~ 0.2 c. 

Since the partonic jet, which produces the Mach wave, 
is created in the very early stage of the collision and 
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FIG. 7: (color online) Mean transverse mass spectra of differ- 
ent particle species (tt", A"*" and protons) at three different 
beam energies compared to data [69l . [73 



traverses the medium until freeze-out, the observed angle 
should also be related to a time average of and not the 
speed of sound at some specific point in time. Because 
the systems spends quite a substantial amount of time 
close to the phase transition region, where the EoS is 
soft, the averaged speed of sound is much lower than the 
limit for an ultrarelativistic gas ^1/3. Still, due to the 
time and space average, it is substantially larger than 
the speed of sound in the transition region. Although 
the experimental result, obtained at much larger beam 
energies (and therefore smaller chemical potentials), is 
not directly related to our results at lower energies, the 
excitation function of (cg) shows a saturation at rather 
moderate energies and, therefore an even lower speed of 
sound in systems created at RHIC seems unlikely. 



RESULTS FOR FINAL PARTICLE PROPERTIES 

In the following, we compare results from our calcu- 
lations with the deconfinement equation of state (solid 
lines) including a first order phase transition, to those 
obtained when an EoS resembling a hadronic resonance 
gas (HG, dotted lines) is applied to the hydrodynamic 
evolution. All results shown are obtained by applying 
the hybrid model to most central (6 < 3.4 fm) heavy 
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FIG. 8: (color online) Excitation functions of the mean trans- 
verse mass of pions (upper panel), protons (upper middle 
panel), negatively charged kaons (lower middle panel) and 
positively charged kaons (lower panel) compared to data 

Iz2,l73,l9l. 



ion reactions (Au+Au/Pb+Pb) in a broad energy range 
from Ei^b = 4 - 160A GeV. 

Note that the observation of particle multiplicity fluc- 
tuations and their kurtosis have become the focus of 
attention concerning the search of the critical endpoint 
[S 0, m , 85 1 . But as has been pointed out in 86] fluc- 
tuations can be very sensitive on the correct treatment 
of conserved quantum numbers like baryon number 
charge and strangeness on an event-by-event basis. It is 
certainly possible, and planned, to analyze these kind of 
event-by-event fluctuations in the current hybrid-model. 
The computational effort for such studies however, is 
much greater than for bulk observables and we therefore 
restricted our present study on non event- by- event 
observables. 



Fig. [5]shows the total yields of different particle species 
as a function of center-of-mass energy of the colliding 
nuclei. One can clearly see that the total multiplici- 
ties of pions and kaons reproduce the experimental data 
(squares) reasonably well. Multistrange hyperons like the 



8 



0.20 
0.15 
0.10 
0.05 
0.00 
-0.05 



-0.10 



■ ■ UrQMD Protons 
■■ UrQMD 71 

— hybrid DE Protons 

— hybrid DE n 

— hybrid HG Protons 

— hybrid HG n 



E, =40 AGeV 

Lab 




♦ Proton Data (NA49) 

* Pion Data (NA 49) 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

y/y^ 




FIG. 9: (color online) The directed flow ui for pions and pro- 
tons at iSiab = 40yl GeV for both hybrid model calculations 
(blue HG and red DE lines) compared to default UrQMD 
results (black dashed lines) and data (green stars and dia- 
monds), yb refers to the beam rapidity [93 ]. 



FIG. 10: (color online) The directed flow vi for protons at 
Eiah ~ 40A GeV for both hybrid model calculations (blue 
triangles HG and red circles DE) compared to default UrQMD 
results (black squares) and data (green stars) . 



E!~ are overestimated at the lowest energies, but follow 
the data very nicely at energies above Ei^h = llA GeV 
(At lower energies strangeness should be suppressed due 
to the chemical non-equilibrium of strangeness). This 
hints to the fact that strangeness is thermalized in heavy 
ion collisions at energies down to i?iab ~ llA GeV. Only 
in the sector of heavy anti-baryons (i.e. S"*") the model 
yields too few particles in the full range of energies inves- 
tigated. One can clearly see that the excitation functions 
of total particle multiplicities do not show any signal of 
the change of the underlying EoS. In contrast, this ob- 
servable seems much more sensitive to the treatment of 
chemical particle freeze out. 

An alternative and more microscopic way to understand 
an enhancement in the production of multi-strange (anti- 
)baryons is related to the fragmentation of the initial 
color flux tubes. While one usually assumes that the ini- 
tially produced color strings fragment independently, a 
high density of strings may allow for a recombination of 
color charges at the string ends, resulting in the forma- 
tion of a fused string with an enhanced color field (known 
as color ropes) 's^, '8l'| . These objects may be seen as an 
alternative to QGP formation. Extensive studies within 
the RQMD model have shown that particle creation from 
these strong color fields results in a substantial enhance- 
ment of multi-strange (anti-)baryons [82,, .83.] . 

Next we turn to the rapidity distributions of different 
particles (pions, kaons and protons). Fig. |6]shows the ra- 
pidity distributions for pions, kaons and protons at three 
different beam energies. For the lowest energy differences 
in the proton rapidity spectra can be observed, although 
both are, within error bars, in agreement with the data. 



For the K'^ and tt^ we observe a good agreement up to 
energies of Slab — 40v4 GeV, while at the highest ener- 
gies (i?iab = 160^ GeV) both hybrid model calculations 
deviate from experiment. This may be related to our 
approximate choice of the hypersurface normal c?cr^^' . 

Since the transverse momentum distributions should 
be sensitive to the transverse dynamics, we next turn 
to the investigation of the distributions of the trans- 
verse mass of different particles (pions, kaons and pro- 
tons). Comparing the hybrid model calculations to data 
(squares), in Fig. [71 we find very good agreement for the 
lower energies. Especially the pion and kaon spectra are 
well in line with the data, while the proton distributions 
deviate slightly. Comparing the momentum distributions 
of the different EoS one again observes no noticeable dif- 
ference between a pure hadron gas EoS and an EoS with 
a deconfinement transition. This hints to the fact that 
the final particle distributions are much more sensitive 
to the freeze-out procedure (thermal distribution of par- 
ticles) and even more to the initial momentum distribu- 
tions given by UrQMD as compared to the effects of the 
expansion dynamics. 

Fig. [51 shows the excitation function of the mean trans- 
verse mass of different particles, at midrapidity. Espe- 
cially for the pions one observes a good agreement with 



[1] Neglecting space like components of the hypersurface leads to 
an overestimate of particle production at large rapidities. As we 
enforce baryon number conservation, the total norm of the rapid- 
ity distribution is conserved and therefore less particles should 
be produced at small rapidities, which can be observed in our 
results for the highest energies 
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data for low energies. At high energies, both hydrody- 
namic calculations overestimate the mean transverse mo- 
menta, while the calculation with the DE EoS gives a 
slightly better description of the data. 

To get an understanding of how much the EoS can 
(ultimately) influence the flnal transverse mass spectra, 
we performed two additional calculations, where we set 
the pressure (as a function of the energy and baryon- 
number densities) to either or p = l/3e throughout the 
whole hydrodynamic evolution. These limits represent 
the stiffest and softest EoS possible. The results from 
this comparison at i^iab = IQOA GeV are depicted in 
Fig. [H as the gray bands, giving an upper (p = l/3e) 
and lower (p = 0) bound on the results from the hybrid 
model. Note that, especially the pion mean transverse 
mass can only be accounted for if the (unphysically) 
soft EoS {p — 0) is applied. Investigations within the 
hybrid model, where an EoS including a flrst order 
phase transition with a very large latent heat (very 
small pressure gradients) is employed, come to a similar 
conclusion 93[. In order to explain transverse flow data 
with ideal fluid dynamics in a hybrid model, an EoS 
with an unphysically large first order phase transition is 
needed. 

The flattening of the mean transverse momentum as a 
function of beam energy was assumed to be a signal for a 
first order phase transition [h^I . But as has been pointed 
out above the slow increase of the mean transverse 
mass at increasing beam energy can be related to non 
equilibrium effects, and may be modeled by the inclusion 
of viscous hydrodynamics. 

Another observable that is deemed to be a signal for 
the deconfinement phase transition is the so called 'anti- 
flow' (or third flow component). In 17, [3 H^l it was 
suggested that a first-order phase transition leads to a 
prominent wiggle in the directed flow vi = (px/pr) as a 
function of rapidity. Data [s^, indeed, show a wiggle in 
the directed flow in most central collisions at Slab = 40^ 
GeV. Fig. [5] shows vi of pions (circles) and protons 
(squares) as a function of rapidity (over yt, the beam ra- 
pidity) for different model calculations, compared to data 
(green stars and diamonds). Here, we compare both hy- 
brid model calculations (with and without a phase transi- 
tion) and a pure transport calculation (UrQMD with- 
out any hydrodynamic stage). As one can see all models 
yield very similar results and none can reproduce the be- 
havior of the proton vi . The pion vi is very well in line 
with data for both hydrodynamic calculations. Note that 
even the sign change for the pions as compared to pro- 
tons is reproduced. Still, changing the EoS does not lead 
to distinguishable differences in the extracted directed 
flow. The directed flow seems to be mostly sensitive to 
the initial conditions and less to the subsequent hydrody- 
namical expansion. This is in line with previous findings, 
where the inclusion of nucleon potentials (mean field ef- 



fects) in transport models improved the description if the 
proton directed flow as measured by experiment 95, 9^ . 



Fig. [To] again shows the directed flow of protons for a 
smaller window in rapidity. The change of sign in proton 
flow can be clearly seen in the data, but neither model 
calculation (hybrid model with and without phase transi- 
tion) describes this phenomena (not even qualitatively). 
Disregarding this single data point, the pure transport 
calculation provides the best description of the directed 
proton flow. 



SUMMARY 



We discussed the influence of chiral restoration and the 
deconfinement phase transition in a hybrid model calcu- 
lation. To model a more realistic equation of state of hot 
and dense nuclear matter we introduced the Polyakov 
loop as the order parameter of deconfinement in a chi- 
ral hadronic model. This EoS was applied to a mi- 
cro-|-macroscopic hybrid model for heavy-ion collisions 
at intermediate energies. Although the phase transition 
leads to a softening of the EoS, we find that the average 
speed of sound in the expanding fireball is of the order 
of Cs ~ 0.3 c over a wide range of energies. This is in 
contrast to findings at the STAR experiment, where 3- 
particle correlation data, in connection with a Mach cone 
scenario, seem to indicate a much smaller average speed 
of sound. 

We compared calculations for final particle properties in 
central A-l-A collisions over a wide energy range to avail- 
able data. Here either a hadron resonance gas EoS or 
the deconfinement EoS was applied. Regarding most ob- 
servables both models give a reasonable description of 
data, without tuning any parameters of the model. Two 
proposed observables for the deconfinement phase tran- 
sition, the flattening of the mean transverse mass ex- 
citation function and the appearance of the third flow 
component (anti-flow), were investigated within our ap- 
proach. We showed that the observed flattening of (ttit) 
cannot be explained by a phase transition, even if a more 
realistic equation of state is applied. Here other mech- 
anisms, like viscosity and other non equilibrium effects 
play an essentially more important role. 
Finally we have shown, that the effect of the EoS on the 
directed flow is almost completely negligible. No anti- 
flow is observed in our calculations and especially pion 
flow is described very well. This leads to the assumption 
that the proton flow is more sensitive to non-equilibrium 
effects. Within this picture, the importance of the (non- 
equilibrium) initial state for flow observables was pointed 
out. 
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